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Abstract. We determined the separations of the radio 
sources in the triangle formed by the BL Lac objects 
1803+784 and 2007+777, and the QSO 1928+738 from in- 
tercontinental interferometric observations carried out in 
November 1991 at the frequencies of 2.3 and 8.4 GHz si- 
multaneously. Difference phase-delay astrometry yielded the 
following separations: 
Aa(i803+784)-(i928+738)=-l''27" 2f 811256±0?000062 
A5(i803+784)-(i928+738)= 4° 30' 2"44833 ±0"00012 
Aa(i928+738)-(2oo7+777)=-0''37"42f503305±0f000033 
A(5(i928+738)-(2007+777)=-3° 54'41"67756 ±0"00013 
Aa(2007+777)-(i803+784)= 2'' 4"45f314561±0f 000052 
A<5(2oo7+777)-(i803+784)=-0° 35'20"77077 ±0^'00013 
We successfully connected differenced phase delays over 7° 
on the sky at 8.4 GHz at an epoch of maximum solar activ- 
ity. The effects of the ionosphere on these VLBI data were 
mostly removed by estimates of the total electron content 
from observations of GPS satellites. The comparison of the 
estimated separation of QSO 1928+738 and BL 2007+777 
with previous such estimates obtained from data at dif- 
ferent epochs leads us to a particular alignment of the 
maps of QSO 1928+738 at different epochs relative to those 
of BL 2007+777, although with significant uncertainty. For 
this alignment, the jet components of QSO 1928+738 show a 
mean proper motion of 0.32±0.10mas/yr and also suggest 
an identification for the position of the core of this radio 
source. 

Key words: astrometry - techniques: interferometric - BL 
Lacertae objects: individual: BL 1803+784, BL 2007+777 - 
quasars: individual: QSO 1928+738 



1. Introduction 

A key need of astronomers is a celestial reference frame. 
Initially they used optical observations to obtain such a 



refere nce frame (see, e.g., Fricke et al. 1988| , Fricke et al 



1991 ). A more rigid and accurate frame can be obtained us- 
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ing Very Long Baseline Interferometry (VLBI) , a technique 
that provides far higher angular resolution (see, for example, 
Johnston et al. 1995). Using differenced phase delays from 



observations of two sources yields relative positions of radio 
sources with submilliarcsecond accuracy and fractional pre- 
cisions reaching ~2 • 10"^ (Guirado et al. [19954 |1995bt Lara 



et al. 1996). Differenced phase delays may provide the most 



precise observable in astrometry. The group-delay observ- 
able, although also capable of providing global astrometry, 
has a statistical standard error larger than the phase-delay 
error by a factor equal to the ratio of the center frequency 
of the VLBI observations to their effective bandwidth. 

In radio astrometric work, it is important to account for 
the propagation medium, including the ionosphere. Sardon 



et al. (1994) showed that the total electron content (TEC) 
of the ionosphere could be determined with high accuracy 
(standard error of about 2%) for all directions from a site 
by using dual-frequency GPS data. The ionospheric contri- 
bution can also be determined from dual-band VLBI ob- 
servations by comparing connected phase delays (or group 
delays) from the two bands. In this paper, we compare those 
approaches for estimating the effect of the ionosphere. 

Our use of a triangle of targets allows an improvement 
with respect to astrometry with pairs of radio sources (e.g., 
Guirado et al. 1995a 1995b, 1998): closure on the sky pro- 



vides a check on the internal consistency of the data reduc- 
tion. 

2. VLBI observations and data analysis 

We observed the radio sources BL 1803+784, 
QSO 1928+738, and BL 2007+777 in right circular polar- 
ization at both 2.3 and 8.4 GHz, on 20 and 21 November 
1991 (MJD 48580 and 48581) from 14:00 to 04:00 UT. 
Using a global VLBI a rray, we recorded with the Mark III 
system (Rogers et al. 1983 ), in modes A (56 MHz) or B 
(28 MHz), depending on the instrumentation at the indi- 
vidual sites. The array consisted of the following antennas 
(in parenthesis, code, diameter, location): Effelsberg (B, 
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100m, Germany); Haystack (K, 37m, Massachusetts, US); 
Medicina (L, 32m, Italy); DSS63 (M, 70m, Spain); and 
phased- VLA (Y, 130 m effective diameter, New Mexico, 
US); Fort Davis (F, 25m, Texas, US); Kitt Peak (T, 25m, 
Arizona, US); Pie Town (P, 25m, New Mexico, US); Los 
Alamos (X, 25m, New Mexico, US); and Green Bank (G, 
43m, West Virginia, US). 

The three radio sources were observed cyclically (cycle 
time of 9 min, clockwise direction) with 2 minutes "on" ra- 
dio source and 1 minute allowed for slewing antennas. All 
sources had total flux densities greater than 1 Jy at the 
epoch of observation. We obtained strong detections from 
each of the baselines for at least 80% of the 2 min "scans" . 

The observations were plagued by problems. Medicina 
and Green Bank had phasing problems. DSS63 recorded 
only 3 hours of data. The phased- VLA recorded only 
8.4 GHz data, and Haystack was set to a wrong sky fre- 
quency. After many fruitless attempts to use part of the 
data from these stations, we decided to use only the data 
from the subset of antennas with very good performance: Ef- 
felsberg. Fort Davis, Pie Town, Kitt Peak, and Los Alamos. 

The data were correlated at the Max-Planck Institut 
fiir Radioastronomie, Bonn, Germany. After fringe-fitting 
the correlator output, we obtained for each observation and 
baseline, estimates of the group delay, the phase-delay rate, 
the fringe phase (modulo 27r), and the correlation ampli- 
tude at each of two selected reference frequencies (2,294.99 
MHz and 8,404.99 MHz). Following standard procedures we 
constructed the visibilities using the calibration information 
provided by staff of the antennas. We made hybrid maps of 
each of the three radio sources at 2.3 and 8.4 GHz using the 
Ca hech package (Pearson |l99l|) and DIFMAP (Shepherd et 
al. 1995 ). Our mapping results are presented in Sect. ^ and 
our astrometric results in Sects. ^ and ||. 

3. Mapping of the radio sources BL 1803+784, 
QSO 1928+738, and BL 2007+777 



A complete sample (S5) (Kiihr et al. 1981 ) of compact ra- 
dio sources has been studied using the VLBI technique over 
the last two decades (Eckart et al. 
al. 



1986, 1987, Witzel et 



These sources were selected from the fifth install- 
ment of the MPIfR-NRAO 5 GHz strong source survey. The 
sources have flat spectra that permit multiband observa- 
tions. All members of this sample have neighbors relatively 
nearby on the sky, allowing phase-reference observations to 



be made (Guirado et al. 1995a, 1998) 



3.1. The BL Lac object 1803+784 

1803+784 is a BL Lacertae type object with z=0.684 and 
V=16.4 (Stickel et al. 1991). It exhibits a radio structure 



with kiloparsec-size components (Antonucci 1986, Strom & 
Biermann |1991| , KoUgaard et al. |1992| , Murphy et al. [1993[ ) 
to the south of the core (at angular distances from the core 
of 2", 37" and 45"). On parsec scales the source sh ows a 
westward j et (E ckart et al. 1986 Pearson & Readhead 1988, 
Schalinski 199C); VLBI monitoring of these milliarcsecond 



(mas)-size components at 22 GHz shows components at an- 
gular distances of about 0.4, 1.4 and 5 mas from the core 



(Eckart et al. |198£, Schahnski 199C, Krichbaum et al. 1990 



1993, 1994a, 1994b). Observations at 22 and 43 GHz also 



reveal other features that appear to be moving between 



the core and the 1.4 mas component. Steffen (1994) has 
modeled this radio source as an adiabatically expanding 
homogeneous plasma jet moving along a helical trajectory. 
The parsec-scale structure is misaligned with the kiloparsec 
structure by '^90°. 
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Fig. 1. Hybrid maps from observations on 20-21 November 
1991 of the BL Lacertae object 1803+784 at 2.3 and 8.4 
GHz. Contours are 0.5, 1, 2, 4, 8, 16, 32, 64, and 90% of the 
peak of brightness, 1.78 Jy/beam, for 2.3 GHz, and -0.25, 
0.25, 0.5, 1, 2, 4, 8, 16, 32, 64, and 90% of the peak of bright- 
ness, 1.40 Jy/beam, for 8.4 GHz. Dotted lines correspond to 
negative contours. The size of the natural restoring beam, 
shown at the bottom left corner of each map, is 2.53x1.97 
mas (PA=58°) for 2.3 GHz and 0.65x0.49 mas (PA=49°) 
for 8.4 GHz. Note that the scales of the two maps are not 
the same. 



In our maps (see Fig. we find well defined components 
to the west of the core, with the jet structure extending up 
to 35 mas from the core at 2.3 GHz and up to 21 mas at 
8.4 GHz. The features in our maps agree well with those 
shown in previous maps. 
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3.2. The QS 01928+738 (4C 73.18) 

QSO 1928+738, also known as 4C 73.18, is a w ell-stu died 
quasar with z=0.302 and ^=15. 5 (Lawrence et al. 1986 ). On 
kiloparsec scales, it exhibits two-sided jet emission feeding 
two symmetrically placed lobes in a north-south orientation 
(Rus k & Rusk |1986| , Johnston et al. |1987| , Murphy et al. 
199^ ). On parsec scales it displays a pronounced one-sided 



southward jet with supcrluminal motion along a sinusoidally 
curved jet (Eckart 1984 Pearson & Readhead 1988, Hummel 
19901 , and Hummel et al. 1992, hereafter p92| ). |H92| describe 



this source as having a helical magnetic field attached to 
a rotating accretion disk. Further, Roos et al. (1993) sug- 
gest that the source consists of two massive black holes, 
with mass ratio smaller than 10, total mass ^lO^M© and 
orbital radius ~10^^cm. Such massive binary black holes 
could be responsible for the sinusoidal jet ridge line which 
has a period of precession of about 2.9 years according to 
the ballistic relativistic jet model. 

Our maps (see Fig. g) show the core-jet structure, with 
components in the jet up to 30 mas south of the core at 
2.3 GHz and up to 13 mas south at 8.4 GHz. 



3.3. The BL Lac object 2007+777 

2007-f 777 is also a BL Lacertae type object with z 
and y=16.5 (Stickel et al 



1991) 



0.342 

It is a core-dominated 
radio source with a westward-directed jet of parsec scale 
with superluminal components discernible up to 6 mas from 
the core (Eckart |198|). This source has also been studied on 
a kiloparsec scale (Antonucci 1986, KoUgaard et al. 1992). 
The kiloparsec structure is aligned with the parsec structure 
and has components '^11" east of the core, and a jetlike 
feature extending 16" to the west with components at ~8.5" 
and ^15.8" from the core. 

For this source our maps (see Fig. ^ show a relatively 
small jet extending 8 mas west of the core. We identify the 
brightest component of the map as the core due to its spec- 
tral power being the same at 2.3 GHz and 8.4 GHz. 



2007+777 



2.3 GHz 



10 - 



- 



-10 - 




10 -10 

Relative R.A. (milliarcsec) 



2007+777 



.4 GHz 



- 



Relative R.A. (milliarcsec) 

Fig. 3. Hybrid maps from observations on 20-21 Novem- 
ber 1991 of the BL Lacertae object 2007-^777 at 2.3 and 
8.4 GHz. Contours are 0.5, 1, 2, 4, 8, 16, 32, 64, and 90% 
of the peak of brightness, 0.69 Jy/beam, for 2.3 GHz, and 
-1, 1, 2, 4, 8, 16, 32, 64, and 90% of the peak of bright- 
ness, 2.06 Jy/beam, for 8.4 GHz. Dotted lines correspond 
to negative contours. The size of the natural restoring beam, 
shown at the bottom left corner of each map, is 2.58x1.97 
mas (PA=87°) for 2.3 GHz and 0.69x0.54 mas (PA=87°) 
for 8.4 GHz. Note that the scales of the two maps are not 
the same. 



4. Astrometric data reduction 

Our astrometric data reduction contains several steps which 
yield nearly "ionosphere-free" and "structure-free" phase- 
connected delays. Using an extensively mod ernize d version 
of software originally written by Robertson (1975), we ana- 
lyze these delays via weighted least-squares to estimate the 
angular separations among the sources. In the next para- 
graphs we discuss in detail each of the steps of this process. 

4-.1. Phase connection 

We used the 8.4 GHz group delays, Tg, together with the 
phase-delay rates, f^, to estimate via weighted- least-squares 
the values of the parameters of our theoretical models in- 
cluding the behavior of the atomic clock at each site rel- 
ative to that at a reference site. Using this astrometric 
model we predicted for each of our ten baselines the num- 
ber of cycles of phase between successive observations of 



the same source to permit us to "connect" the phase de- 



lays (Shapiro et al. 1979). Denoting the reference frequency 
(8,404.99 MHz) by Vr and the cycle time of the observations 
(9 min) by Ai, we can represent the constraint on the stan- 
dard error, cr{T^), needed to insure proper phase connection: 
c(f0) <C (t'rAt)"^, i.e., cr(T0) must be small compared to 
0.2ps-s-\ In our observations, this criterion is satisfied at 
8.4 GHz for our observations and essentially complete phase 
connection was feasible. However, phase connection was not 
possible for the data at 2.3 GHz, due to the strong effect of 
the ionosphere and the consequent high scatter in the phase- 
delay data, leading to a violation of the criterion stated 
above. 

We also considered possible overall offsets (equivalent to 
an integral number of cycles of phase for all the baselines 
involving a given site and source, i.e., constant clock off- 
sets or fringe ambiguity offsets). To estimate such possible 
offsets, we considered the fringe ambiguity offsets as param- 



E. Ros et al.: High precision astrometry on B1803+784/Q1928+738/B2007+777 



1 1 1 1 1 1 
1928+738 


1 1 1 1 1 
2.3 GHz 




m - 






will 




c 


) : 

O - 








1928+738 



8.4 GHz 




O 



10 -10 

Relative R.A. (milllarcsec) 



5 0-5 
Relative R.A. (milllarcsec) 



Fig. 2. Hybrid maps from obser- 
vations on 20-21 November 1991 
of the QS01928-f738 (4C 73.18) 
at 2.3 and 8.4 GHz. Contours are 
0.5, 1, 2, 4, 8, 16, 32, 64, and 
90% of the peak of brightness, 
1.16 Jy/beam, for 2.3 GHz, and 
-0.5, 0.5, 1, 2, 4, 8, 16, 32, 64, 
and 90% of the peak of bright- 
ness, 1.18 Jy/beam, for 8.4 GHz. 
Dotted hnes correspond to neg- 
ative contours. The size of the 
natural restoring beam, shown at 
the bottom left corner of each 
map, is 1.98x2.62 mas (PA=80°) 
for 2.3 GHz and 0.72x0.53 mas 
(PA=76°) for 8.4 GHz. Note that 
the scales of the two maps are not 
the same. 



eters whose values were to be estimated and we introduced 
one of them into the global weighted-least-squares analysis 
for each site for each radio source. The values obtained were 
all less than 24 ps, a fifth of an ambiguity interval. We con- 
cluded that the fringe ambiguity offsets for the phase-delay 
data at 8.4 GHz are all zero. 

4-. 2. Radio source structure correction 

The selection of the reference point in each map is a criti- 
cal task. The meaningful accuracy of our final sky-position 
estimate can be no better than the accuracy with which we 
can identify a specific "physical" reference point in our map. 



Using the Caltech package (Pearson 1991), we made special 
renditions of the images of our radio sources with a pixel size 
of 0.01 mas for 8.4 GHz (about 70 times smaller than the size 
of the CLEAN beam). For Figs. |l|-^, the delta- function com- 
ponents from the CLEAN procedure had been convolved 
with the corresponding Gaussian beams. For these special 
renditions of the images we also convolved the CLEAN com- 
ponents with their corresponding Gaussian beams. We chose 
as the reference point the center of the pixel corresponding 
to the peak of the brightness distribution, for each radio 
source and frequency. For a source with a "sharp" peak, 
the standard error in this location is proportional to the 
ratio between the size of the interferometric beam and the 
SNR of the peak of brightness in the map (see, for example, 
Fomalont 1989). Using our reference points, we calculated 
the structure contributions to the phase delay. None of the 
contributions is larger than 25 ps. Thus, the structure con- 
tributions at 8.4 GHz are ~0. 2 cycles (no more than 70°) for 
QSO 1928-1-738, the source with the largest structure. 



with flux densities larger than 25% of the component with 
the highest flux density, in the area covered by the "main 
lobe" of the synthesized beam. The difference in the posi- 
tions determined from these two criteria is < 20 fias. 

4-3. Troposphere correction 

The neutral atmosphere introduces a delay Tatm in the op- 
tical path length of the radio wave of about 2.3 m in the 
zenith direction. From meteorological data obtained at each 
antenna site, we calculated (see, e.g., Saastamoinen 1973| ) 
a priori values for the two components of the delay at lo- 
cal zenith: the dry and the wet. We show these values in 
Table I. 

For the weighted-least-squares analysis, we represented 
the zenith delay as a continuous, piecewise-linear function of 
time and determined it from meteorological values sampled 
every two hours during the observation. We used the Chao 



Another choice (see Guirado et al. 1995a) of reference 
point is the centroid of those delta-function components 



mapping function (Chao 1972, 1974) to approximate the 
delay at elevation angles other than the zenith. These angles 
were always greater than 20° , and for these delays the Chao 
mapping function yields results consistent with those from 
other more accurate mapping functions. 

4.4- Ionosphere correction 

The plasma delay due to the Earth's ionosphere was the 
most difficult contribution to estimate. With observations 
made at two frequency bands simultaneously, the iono- 
spheric contribution can largely be removed from some 
observations by taking advantage of the dependence 
of the effect. The delay of a radio wave propagating 
through a (weak) plasma can be approximated by: A rjon — 
±{kEc) / [cv^) (see, for example, Thompson et al. 1986), 
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Table 1. Parameters of the theoretical model used in the weighted-least-squares analysis. 



Source coordinates (J2000.0)'' 

1928+738'' a = 19''27™48f 495213 5 = 73°58'01"56995 

Antenna site coordinates and propagation medium parameters'' 



Radio 








Axis 


Mean tropospheric 


telescopes'' 


Cartesian coordinates 


[m] 


offset" [m] 


zenith delay' [ns] 




X 


Y 


Z 




Tdry 


"^wct 


B 


4033947.528 


486990.428 


4900430.727 


0.0 


7.4 


0.2 


F 


-1324009.072 


-5332181.961 


3231962.515 


2.123 


6.4 


0.4 


P 


-1640953.650 


-5014816.009 


3575411.928 


2.136 


5.9 


0.1 


T 


-1995678.575 


-5037317.707 


3357328.183 


2.138 


6.2 


0.0 


X 


-1449752.303 


-4975298.589 


3709123.975 


2.139 


6.2 


0.0 



Earth Tides 

Radial Love number, h=0. 60967 

Earth Orientation Parameters" 
Nutation 

lAU 1980 Nutation Series + 
lERS corrections for MJD=48580: 
dV = -0!'01238 de = -0"00201 



Horizontal Love number, 1=0.085 Tidal lag angle, 9 — 0?0 



Polar Motion and UTl 

Wobble: x=0"25221 y=0''27880 
lERS LOD correction= 0.002422s 
UT1-UTC= -0.027891s 



Precession Constant (J2000.0)«5 

p=5029".0966/Julian century 
Mean obliquity (J2000.0)'5 

eo=23°26'2l"448=8438l"448 



^ The reference coordinates are taken from 1993 lERS Annual Report (lERS 1994); they are used only for convenience. We used 
standard deviations from the lERS parameter estimates, but increased them twofold to be conservative. These increased standard 
deviations are: 0.035 ms and 0.28 mas for the {a, 5) QSO 1928+738 coordinates, 2 cm in each coordinate of site position, Imas in 
each pole coordinate, and 0.04 ms in UTl-UTC. See Table |. 
The position of QSO 1928+738 was held fixed and the positions of the other two sources were determined with respect to it. 

'^ The site positions were calculated by linear extrapolation (using the constant velocities for the coordinates obtained from lERS 
(1994) and applied to its 1993.0 positions). 

'' See text for key. 

" All the antennas are altitude-azimuth. 

* Mean values of th e zeni th delays, computed using the meteorological values provided by the staff of the observing antennas (see, 
Saastamo inen 1973 ) . 



e.g. 

= Lieske et al. 



(1977) 



where k—AO.3 m'^s^^, c is the speed of light (m-s^^), ly 
the frequency (Hz), and Ec the integrated electron content 
(m^^) along the line of sight; the plus and minus signs re- 
fer to group and phase delays, respectively. Unfortunately, 
for our observations, which were made near a maximum 
in solar activity, we could not phase connect the delays at 
2.3 GHz, as noted in Sect. 4.1. Although the group-delay 
data at 2.3 GHz and at 8.4 GHz can be combined to yield 
ionospheric corrections, these are not accurate enough to 
allow us to phase connect the 2.3 GHz delay data. 

Because we believed the corrections for the ionosphere 
would be somewhat more accurate, we used estimates of 
Ec from Global Positioning System (GPS) data (Sardon 
et al. 1994) instead of from the dual-band VLBI data. 



The GPS system uses two frequency bands, one each at 
l,575.42MHz and l,227.60MHz. From data obtained si- 
multaneously at one site from different satellites, we esti- 
mated the total electron content in the zenith direction. 
With such estimates from a number of sites, together with 
our model of the electron distribution in the ionosphere, 
we can interpolate to obtain the delay in a given direc- 
tion for a given (relatively nearby) site (see Klobuchar 
1975). We therefore obtained data from several GPS sta- 



(CA, US, 116.890''W, 35.426°N) and Pinyon Flats (CA, US, 
116.458°W, 33.612°N) available at the GARNER archives 
at SOPAC, University of California-San Diego, to estimate 
the TEC for each of our VLBI observations for the follow- 
ing antennas: Fort Davis (103.944°W, 30.635''N), Pie Town 
(108.119°W, 34.301°N), Kitt Peak (in.612''W, 31.956°N), 
and Los Alamos (106.245''W, 35.775''N), in each case in- 
troducing a longitude (time) correction and an appropri- 
ate "slant" factor. In our relatively crude procedure, fol- 
lowing the geometry from Klobuchar (1975), we assumed 
that the TEC was concentrated at an altitude of 350 km. To 
determine the TEC for Effelsberg (6.844°E, 50.336°N), we 
used the data from Herstmonceux (UK, 0.336''E, 50.867°N) 
and WettzeU (Koctzling, Germany, 12.789°E, 49.144''N) ob- 
tained from the archives mentioned above. Further details 



tions for the time span of our observations. For the U.S. 
part of the array we used the TEC data from Goldstone 



will be provided in Ros et al. (1999). We used the GPS re- 
sults as described above to deduce ionosphere corrections 
for the line of sight from each VLBI station to each ra- 
dio source for each scan, and to thereby correct the group 
delays (at both frequencies) and phase delays (at 8.4 GHz). 
The corrections were as large as 10 turns (cycles) at 8.4 GHz 
for some scans on intercontinental baselines. The statistical 
standard errors of the GPS corrections at 8.4 GHz are 30 ps, 
a fourth of an ambiguity interval, corresponding to a TEC 
standard error of 1.5x 10^^ el/m^. 
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4-5. Differenced observables 

After the phase-connected observables at 8.4 GHz for each 
source were corrected for source structure and ionospheric 
effects, we formed the differenced phase delays by subtract- 
ing the phase delay of each observation of a radio source 
of the triangle from the previous one observed in the cycle. 
The differenced phase delays are thus partially free from 
unmodeled effects that remain in the phase delays for each 
one of the sources, such as those effects due to instabilities 
of the time standards used at the sites and inadequacies 
in the modeling of the troposphere. When working with 
sources nearby to one another on the sky, subtraction of 
results from consecutive observations tends to cancel errors 
since the received signals from the two sources follow simi- 
lar paths at nearly the same epochs. With sources at larger 
angular separations, such tendencies to cancel are weaker: 
not only are the spatial separations of the ray paths larger, 
but the epochs of successive observations tend to be further 
apart, due to the limited slew speeds of the antennas. In our 
case the source separations range from 4° to 7° . 

Next, the undifferenced and differenced phase-delay data 
from the three radio sources were used in a weighted-least- 
squares analysis (the use of different subsets of data, e.g., 
undifferenced data from only one or two of the three radio 
sources, or the use of two or three sets of differenced data 
does not alter the results or the statistical standard errors 
significantly). We included undifferenced observables to bet- 
ter estimate the long-term behavior of the station clocks. For 
each baseline, we scaled separately the standard deviations 
of the undifferenced phase delays for each radio source and 
the differenced phase delays for each pair of radio sources 
to give a root-weighted-mean-square of unity for the corre- 
sponding set of postfit residuals. The resulting undifferenced 
phase delays have about threefold larger standard errors 
than do the differenced phase delays. 

At this stage of the data reduction we again checked for 
possible nonzero fringe ambiguity offsets and possible errors 
in phase connection. Estimating the fringe ambiguity offsets 
in a least-squares analysis again yielded small values (a few 
picoseconds), showing that no changes were necessary. 

5. Angular separations 

With these phase-connected, nearly ionosphere- and struc- 
ture-free phase-delay data, we estimated the angular sep- 
arations among the pairs of radio sources of the triangle. 
Specifically, we estimated the coordinates of BL 1803-1-784 
and BL 2007-f 777 with respect to the fixed coordinates of 
QSO 1928+738. We used the data of the three radio sources 
jointly to obtain these solutions. The results arc shown on 
Table ^ and the postfit residuals of the differenced phase 
delays in Fig. ^. 

We did not estimate in the weighted-least-squares anal- 
yses the tropospheric zenith delays, the site coordinates, the 
coordinates of the reference source (i.e., QSO 1928-1-738; us- 
ing either of the other two sources as reference does not alter 
the results significantly), the coordinates of the Earth's pole. 



Table 2. Relative right ascensions and declinations in 
J2000.0 coordinates of the radio sources BL 1803-1-784, 
QSO 1928-^738 and BL 2007-^777, labeled as A, B, and C, 
respectively, obtained by a weighted-least-squares analysis 
of the connected phase delays. 



Angular separations'* 


A- 


-B 


Aq 
AS 


-l'^27'" 2f811256±0f000062 
4° 30' 2!'44833 ±0!'00012 


B- 


-C 


Aq 
A<5 


-0" 37'"42f 503305±0f 000033 
-3° 54'41^'67756 ±0"00013 


C- 


-A 


Aq 
AS 


2" 4'"45f314561±0f000052 
-0° 35'20!'77077 ±0!'00013 



^ Estimated errors are based on the combination of the various 
standard errors included in the analysis: statistical standard 
errors, and estimated standard errors in the values of the 
(fixed) parameters of the theoretical model, in the ionospheric 
model from TEC measurements, and in the reference point 
chosen for each map. See text and Tables ^ and ^. Correlations 
among the estimated coordinates for sources A and C, and 
among the sensitivities of those estimates to the parameters 
listed in Table ^, have been taken into account in computing 
the standard errors of the C-A separation. 



and UTl-UTC. Instead, to accommodate software limita- 
tions we fixed their values equal to their a priori values (see 
Table ^ , and accounted for their a priori standard errors in 
our subsequent sensitivity analysis (see Table ^). For each 
solution of the angular separations we estimated the overall 
standard errors of the components of the relative position 
between the radio sources as the root-sum-square (rss) of 
the standard deviations of all the contributions discussed 
above (see Tables || and ^ . 

We define "sky closure" as the sums of the compo- 
nents of the separation vectors between the three radio 
sources determined pairwise. For right ascension a, we 

have: Ca = J2^o^ = ^"(c - b) + ^^^b - A) + ^a(A - c), 
where A, B, and C stand, respectively, for BL 1803-1-784, 
QSO 1928-f738, and BL 2007-1-777. Closure in declination 5 
is defined correspondingly. The results of the sky-closure test 
are Ca = —66 ± 360 /ias and = 15 ± 220 /ias, that is, zero 
to well within the indicated standard errors. We estimated 
these errors as the root-sum-square of the standard errors 
of the determinations of the pairwise separations. Conse- 
quently, these sky closure errors are upper bounds since we 
ignored the correlations among the pairwise estimates of 
separation. 

6. The pair QSO 1928+738/BL 2007+777: 
registrations and proper motions 

We have obtained the pairwise separations among the radio 
sources BL 1803+784, QSO 1928+738, and BL 2007+777 
from observations taken on 20-21 November 1991 through 
a weighted-least-squares analysis of undifferenced and dif- 
ferenced, connected, nearly ionosphere- and structure-free 
phase-delays, r^. 



7 



1928+738 - 1803+784 2007+777 - 1928+738 1803+784 - 2007+777 



100 r 

50 

-50 

-100 i_ 
100 

50 - 

- 
-50 

-100 i_ 
100 
50 

-50 

-100 i_ 
100 
50 

- 

" -50 - 
^ -100 -_ 

^ 100- 
^ 50 - 

oj - 

J ' 

-100 -_ 
■g 100- 
g 50- 

s 0- 

S -50 - 

2 -100 -_ 

° 100 - 

2 50 - 

-g - 

I -50 - 

I -100- 
I 100 
£ 50 

-50 

-100 i_ 
100 
50 

-50 

-100 i_ 
100 
50 

-50 

-100 L 



15 



BF 



BF 



BP 



BT 



BX 



BP 



BT 



BX 



FP 



FP 



111 ,1 Ji 



FT 



FT 



FX 



FT 



XP 



FX 



FT 



XP 



XT 



XT 



BF 



BP 



BT 



BX 



FP 



FT 



FX 



PT 



XP 



XT 



20 



25 



15 



20 25 
UT (hours) 



15 



20 



25 



Fig. 4. Postfit residuals of the 

differenced phase delays for all 
baselines and source pairs in- 
volved in our analysis (see text). 
Full (vertical) scale is ± one 
phase cycle (±27r equivalent), 
i.e., ±119ps. 



Table 4. Error contributions to the three radio source positions obtained from phase delays. 



1803+784 1928+738 2007+777 





6Aa 


SA5 


5Aa 


5A6 


SAa 


SA5 


Cause 


[Mas]" 


[/ias] 


[//as]* 


[/uas] 


[/ias]" 


[/xas] 


Statistical standard error 


12 


11 


_b 


_b 


12 


9 


Fixed parameters in astrometric modeF 


179 


118 


_b 


_b 


98 


123 


Reference point in the map'' 


14 


14 


23 


12 


23 


12 


TEC 


42 


21 


_b 


_b 


16 


14 


rss (estimated standard error) 


186 


121 


23 


12 


103 


125 



" Wc use the factor 15- cos 5 to convert /ts to /ias, wliere 5 is tlie declination of the radio source. 

^ We toolc QSO 1928+738 as reference and estimated the positions of the other two sources relative to the position of this reference. 
" See Table 3. 

Each tabulated value includes the estimated standard error in determining the peak of brightness due to the noise in the map 
and an additional contribution equal in magnitude to the differences between the location of this peak and that of the centroid of 
the brightest pixel (see text). 

" The tabulated values are deduced from a sensitivity analysis for an estimated TEC standard error of 1.5x10^® el/m^ at each site; 
we assume that the TEC errors for each site are uncorrelated with those for the other sites. 



8 



E. Ros et al.: High precision astrometry on B1803+784/Q1928+738/B2007+777 



Table 3. Contributions SAa and SAS to the standard errors 
of the estimates of the sky coordinates of QSO 1928+738 rel- 
ative to those of each of the two other radio sources obtained 
from the error sensitivity study. 

Sensitivity analysis results^ 

1803+784 2007+777 
Standard SAa SAS SAa SAS 



Parameter 




deviation 


[Mas]"^ 


[/ias] 




[/ias] 


Coordinates 


a 


145 ^as'^ 


30 


53 


35 


24 


of 1928+738"^ 


S 


280 fias 


96 


20 


47 


12 


Tropospheric 


B 


0.1ns 


105 


20 


22 


105 


delay 


F 


0.1ns 


48 


27 


22 


14 


at the zenith 


P 


0.1ns 


30 


22 


22 


6 


of the radio 


T 


0.1ns 


12 


9 


22 


7 


telescopes 


X 


0.1ns 


27 


23 


19 


7 




B 


2 cm 


60 


62 


41 


40 


Radio 


F 


2 cm 


27 


18 


16 


16 


telescope 


P 


2 cm 


9 


16 


13 


12 


coordinates'^ 


T 


2 cm 


12 


18 


9 


10 




X 


2 cm 


33 


15 


13 


13 


Earth 


X 


1 mas 


3 


2 





1 


pole'' 


y 


1 mas 











1 


UTl-UTC'' 




0.04 ms 


33 


60 


38 


27 


rss'' (estimated 














standard error) 






179 


118 


98 


123 



Each entry in the last four columns was obtained by altering 
the a priori value of the parameter in turn by one standard 
deviation and repeating the least-squares analysis with the 
resultant table entry being the absolute value of the difference 
between each result and the "nominal" one. 
^ We use the factor 15- cos (5 to convert /is to /ias, where S is the 
declination of the radio source. 

Stan dard deviations are twice the values published in lERS 
( [1994| ) (see Table §). 

The root-sum-square (rss) of the values in the corresponding 
columns. 



Considering previous observations in addition to ours, 
we investigated apparent changes in the relative posi- 
tions of brightness features in the sources QSO 1928+738 
and BL 2007+777. These radio sources were previously ob- 



served astrometrically in 1985.77 (Guirado et al. 1995a) and 
1988.83 (Guirado et al. |1998[ ). For them, we can compare the 
estimated separations for the three epochs. This comparison 
is shown in Table ^ and Fig. |[ For each epoch, the refer- 
ence position for each source is the peak of brightness of 
the corresponding map. Our result from epoch 1991.89 lies 
between the results from the epochs 1985.77 and 1988.83. 

We can obtain some insight into this result by investigat- 
ing the alignment of the maps. The jets of QSO 1928+738 
and BL 2007+777 are oriented nearly orthogonally on the 
sky: The jet of QSO 1928+738 is directed to the south and 
that of BL 2007+777 to the west. This near orthogonality of 
the structures was one of the reasons we selected this pair. 
Indeed, if the cores of the radio sources are assumed station- 
ary as in the standard model, the motions in right ascension 
and declination of the jet components can be studied with- 



es 

B 
to 

< 



0.8 


: 1928+738/2007+777 


1 1 1 1 1 1 1 1 1 ^ 
- 


0.6 






SX'88 


0.4 




X'91 


- 


0.2 








0.0 




C'85 




0.2 








0.4 






< < . 1 . . . 1 . ~ 



0.6 0.4 0.2 0.0 -0.2 -0.4 -0.6 
Aa cos5 (mas) 

Fig. 5. The estimates of the coordinates of QSO 1928+738 
minus those of BL 2007+777 (see Table 5). The origin is 
chosen as in Table 5. Points are labeled to indicate the fre- 
quency band(s) used for phase-delay astrometry and the 
year of each epoch. S corresponds to 2.3 GHz, C to 5 GHz, 
and X to 8.4 GHz. SX indicates that both S and X bands 
were used simultaneously. 



out confusion, using only the differenced astrometric results. 
Alternatively, any apparent motion between the cores that 
cannot be plausibly assigned to some other cause must bring 
the standard model into question and/or the identification 
of the core. For example, the possible emission of a new, un- 
resolved component along a jet direction could change the 
position of the peak of brightness of the corresponding radio 
source. 

An ejection of a component in QSO 1928+738 might im- 
ply a declination-oriented shift of the peak of brightness and 
might affect the difference in declination between this ra- 
dio source and BL 2007+777, but it would be unlikely to 
affect significantly the difference in right ascension. A sim- 
ilar analysis holds for BL 2007+777. Therefore, we inter- 
pret the predominantly north-south changes (see Table |^) in 
the separation of QSO 1928+738 and BL 2007+777 between 
these epochs as shifts with respect to the chosen reference 
position which are caused by changes in the structure of 
QSO 1928+738. In other words, the chosen reference posi- 
tion is not the (center of the) true core. In this scenario, the 
maps of QSO 1928+738 are responsible for all changes: shifts 
of -170, 530, and 260 /ias in declination, and 20, —190, and 
30 /ias in right ascension for epochs 1985.77, 1988.83, and 
1991.89, respectively. Accordingly, the corresponding align- 
ment of the maps of BL 2007+777 (Fig. ^ shows no change 
in the position of the chosen reference point. The maps have 
been overresolved to show more sharply the reference points 
chosen, which have been placed at the origins of the maps. 
We show in Fig. the steps followed to align the maps of 
QSO 1928+738 and the resulting registration. On the top. 
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Table 5. Estimates at each of three epochs of the J2000.0 coordinates, Aa and AS, of QSO 1928+738 minus those of 



BL 20074 


-777^^. 








Epoch 


Antennas'' 


Frequency 
[GHz] 


Aa + 37'"42f 503330'' 
[ms] [mas]'^ 


A(5 + 3°54'41"67780'' 

[mas] 


1985.77 
1988.83 
1991.89 


B,L,F,K,0 
B,L,F,M,t 
B,F,P,T,X 


5 

8.4&2.3 
8.4 


0.005±0.053 0.02±0.22 
-0.045±0.042 -0.19±0.17 
0.008±0.025 0.03±0.10 


-0.17±0.15 
0.53±0.13 
0.26±0.13 



^ Since we assign the shifts in the separation between the two radio sources to QSO 1928+738, we present here the results obtained 
estima ting t he position of QSO 1928+738, while that of BL 2007+777 is h eld fix ed at a = 20''05"'30f 998546, S = 77°52'43!' 24777 
(lERS 1994). The reference coordinates are the same as in Guirado et al. (199^). 

^ The symbols correspond to the following antennas (with diameter and location given in parentheses): B, Effelsberg (100 m, 
Germany); L, Medicina (32m, Italy); K, Haystack (37m, Massachusetts); O, Owens Valley (40m, California); M, DSS63 (70m, 
Spain); F, Fort Davis (25m, Texas); P, Pie Town (25m, New Mexico); t, Onsala (20m, Sweden); T, Kitt Peak (25m, Arizona); X, 
Los Alamos (25m, New Mexico). 

To facilitate comparisons between the values of Aa and AS in the table, we used the factor 15 • cos 6, where 6 is the declination of 
QSO 1928+738, to convert ms to mas. 



we have the three maps aligned at the peaks of brightness 
as originally determined. In the middle we show the result 
of shifting those maps based on the results from astrome- 
try and the assumption of no core motion in BL 2007+777. 
At the bottom, we show the same alignment but with the 
images overresolved to sharpen the distinction between the 
components, albeit at the risk of being misled by the over- 
resolution. 

The remarkable lack of flux density from the northern- 
most part of the maps from the first and third epochs, com- 
pared to that from the second is hard to understand and 
represents a clear puzzle in this scenario. Further, it is dif- 
ficult to understand the proposed shift, with respect to the 
relative-position reference shown in Table ^, in the align- 
ment (Fig. I) of the maps of QSO 1928+738 from 1985.77 
and 1988.83 in terms of a motion of particular features in 
the map (Guirado et al. 1998| ), another weakness in this 



scenario. Nevertheless, following through, we compare these 
two maps with that from 1991.89 and conclude tentatively 
that at each of the three epochs the core is centered about 
1.5-2 mas north of the previously chosen reference point. 
This new identification of the core corresponds to the north- 
ernmost feature of the 1988.83 map (undetected in this sce- 
nario in 1985.77 and 1991.89). 

The core of QSO 1928+738 seems to be continually eject- 
ing components, at an estimated mean rate of one every 1.6 
years ( H92| ). The shifts displayed in Fig. |^ are correspond- 
ingly "explained" in terms of motions of the peak of bright- 
ness, due to newly emerging components. In this scenario, 
we assume that our originally selected reference point cor- 
responded to a different jet component in each map. The 
changes in the flux distribution over six years, a period in 



which, according to EI92, up to four components are likely 



to have emerged, might well have caused the peak of bright- 
ness to have moved during this period. For the 1985.77 map 
at 5 GHz, a difference in the location of the peak of bright- 
ness and that of the core of the radio source could also be 
due in part to opacity effects, but with a value for the shift 
likely to be less than 0.2 mas (applying a dependence with 



kX>^ (Marcaide et al. |1985| ), with 0.7 < /? < 2, and k cor- 
responding to 0.7 mas for A3. 6-13 cm, we obtain 0.10 mas 
(/3 = 2) to 0.21 mas (/3 = 0.7) for A3.6-6cm). This upper 
limit is small compared with the apparent astrometric shift. 

A comparison of this Fig. |^ alignment with the results 
provided by H92 provides new insights into the problem. 
H9| identifled the core of QSO 1928+738 as the north- 



ernmost (and most compact) component of the maps at 
22 GHz, based on a 5 year monitoring program. These au- 
thors reported that this core had a flat sp ectru m, because of 
the similar 8.4 and 22 GHz results. Thus H92 aligned their 
maps with this core to determine the proper motions of the 
jet components that they reported. They also developed a 
model for the kinematics of the jet components involving su- 
pcrluminal motion along a southward traveling sinusoidally 
curved jet ridge line. In the model this path results from 
a 2.9 year period of precession of a binary black hole sys- 
tem that would produce a sine function in the maps with 
a "wavelength" of about 1.06 mas and an "amplitude" of 
about 0.09 mas. This predicted amplitude (in right ascen- 
sion) is of the same order of magnitude as the standard 
error of our astrometric results. The set of components in 
the II92 model defines a "traveling sine wave" moving at 



0.28+0.05 mas/yr, with no significant motions of the com- 
ponents with respect to that of the sine wave. 

Our observations do not reach the resolution obtained by 



H92 , due primarily to our almost threefold lower radio fre- 



quency. We therefore cannot discern the inner structure of 
the radio sources in as much detail as EI92| , but we have been 
able to follow the changes in location of the brightest fea- 
tures by means of our astrometry. By overresolving our map 
and the maps of Guirado et al. ( |1998| ), we tried to identify 
the components previously reported by prI92 . Specifically, we 
inferred that four jet components were emitted between our 
first and third sets of observations, consistent with the re- 

These four components are 



suits and predictions of H92 



labeled from south to north in Fig. 
A3. 



as AO, Al, A2, and 
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Fig. 7. Illustration of steps fol- 
lowed to try to properly align the 
maps. The astrometric shifts in 
milliarcseconds are represented 
(not to scale) at the bottom 
of the middle figures. The hor- 
izontal dotted line in the mid- 
dle figures passes through the 
new map origin. The astrometric 
shifts are the displacements from 
the origins originally chosen for 
the maps (upper figures) to the 
new origins derived from astrom- 
etry (see Table 5). 



AO is identified tentatively as the oldest of these new 
components (this component was labeled B by H92 , as indi- 
cated in the upper part of Fig. ^ . AO can also be identified 
in the maps of f 988.83 and f 991. 89 in Fig. |. It is fainter 
in these two maps than in the one of 1985.77. Our estimate 
of its proper motion, based on our conjectured astrometric 
alignment, is 0.34±0.05mas/yr, with the quoted error being 
the statistical standard error (see Table|6|). 

Al is a component that emerged between 1985.77 and 
1988.83. It corresponds to the peak of brightness in the map 
of 1988.83 and can also be identified in the 1991.89 map, 



yielding an estimated proper motion of 0.32±0.05mas/yr 
(see Table |). 

A2 emerged later and represents the peak of brightness 
in the 1991.89 map; in 1988.83 it had been near our initially 
assumed position of the core. Our estimate of its proper 
motion is 0.27±0.05mas/yr (see Table||). 

The newest component, A3, apparently emerged after 
the 1988.83 epoch. It is, we believe, south of the position 
of the (undetected) core, but, of course, north of the rest of 
the structure. 



Component C marked in the EI92 maps can be identified 
in our maps at 5 and 8.4 GHz. The different relative bright- 
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22 GHz 1988.73 1989.71 

1985.10 1986.40 1987.44 1989.25 




1988.83 1991.89 
8.4 GHz 8.4 GHz 

Fig. 8. Comparison, on the same angular size and temporal scales, of our results and those of H92 on QSO 1928+738. 
At the top, maps from H92's 22 GHz observations restored with circular restoring beams of 0.25, 0.30, 0.30, 0.25, 0.25, 
and 0.25 mas diameter, respectively (left to right). At the bottom, maps from our alignment, with restoring beams of 
0.3, 0.3, and 0.4 mas diameter (left to right). We propose identifications of corresponding features and components: For 
example, H92's C can be identified at 1985.77 and 1991.89. The maps at 1988.73 (this paper) and 1988.83 (H92) show 
similar features. We also propose identifications of components Al and A2 in both sets of maps. The dashed and solid 
lines indicate the proper motions of the components. The horizontal (dotted) lines are drawn for reference. 



1985.77 
5 GHz 



Table 6. Position*^ and proper motions of components in 
the maps of QSO 1928+738. 







Epoch 








1985.77 


1988.83 


1991.89 


Proper 




5 GHz 


8.4 GHz 


8.4 GHz 


motion 




[mas] 


[mas] 


[mas] 


[mas/yr] 


A2 




-1.09+0.10 


-0.25+0.10 


0.27+0.05 


Al 




-0.53+0.10 


0.45+0.15 


0.32+0.05 


AO 


-1.12+0.30 


-0.06+0.15 


0.98+0.15 


0.34+0.05 


C 


-0.49+0.15 


0.44+0.10 


1.43+0.15 


0.32+0.03 


D 


0.08+0.15 




2.31+0.15 


0.36+0.04 



* Angular distances of components are given with respect to the 
origin chosen, based on the astrometric alignment of the maps 
(see text and Fig. 7). 



nesses of components C and D in the 22 and 5 GHz maps of 
1985 might be due to C having still been self-absorbed at 
5 GHz while D was already in the optically thin phase of its 
evolution. Our estimate of the proper motion of component 
C from the three maps is 0.32+0.03 mas/yr (sec Table 

It is remarkable that D is not visible in the 1988.83 map, 
even though it is the brightest feature of the 1985.77 image. 



Similarly, C is at least twofold brighter and also more ex- 
tended in the 1988.83 map than we would have expected 
from its evolution along a jet. The brightness change and 
the extension change can perhaps be explained by the un- 
dersampling of the (w, v) plane for intermediate spatial fre- 
quencies. Table I shows that the array at this epoch con- 
sisted of four European radio telescopes and one in the U.S. 
(Fort Davis) . This array provided sampling only at 30-50 MA 
for six baselines and at about 200 MA for the four inter- 
continental ones. Intermediate spatial frequencies were not 
sampled. Hybrid mapping could have assigned to the inner 
features the flux density of these undetected components, as 
model fitting confirms. Specifically, for this poor sampling 
of those intermediate spatial frequencies, the visibilities of 
the map from Guirado et al. (1998) are essentially equiva- 
lent to those generated by a model, consisting of elhptical 
Gaussian components fit to the observed visibility data, that 
includes a component located where we believe component 
D is located. The other components of the model were lo- 
cated as found in the hybrid maps. The proper motion of 
component D, deduced from only our first and third epochs, 
is 0.36+0.04 mas/yr (statistical standard error). 

Our proper motion results agree with those reported by 



H92. Although we have lower resolution data than they. 
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-2 -4 -6 

Relative R.A. (milliorcsec) 

Fig. 6. Alignment of the maps of BL 2007+777 based on 
the astrometric results. Maps are overresolved and contours 
represent from 2 to 90% of the peak of brightness. Con- 
volving beams (shown at left on each map) are circular, 
with diameters of 0.75, 0.5, and 0.5 mas, top to bottom. 
The crossings of the dashed lines show the reference points 
used in the astrometry. The apparent differences in sepa- 
ration at different epochs between the two radio sources 
QSO 1928-1-783 and BL2007-f777 have been entirely at- 
tributed to QSO 1928-f 738 in this scenario (see text). 



we took advantage of our ability to align the maps of 
QSO 1928-1-738 based on our accurate astrometric results, 
to obtain a longer time base. 

7. Conclusions 

We have shown that at 8.4 GHz phase delays can be "con- 
nected" reliably for sources 7° apart on the sky even at 



epochs of maximum solar activity. Guirado et al. (1995a) 
and Lara et al. (1996) demonstrated such phase connection 
for angular distances of ~5°, but not for epochs of high solar 
activity and rapidly varying TEC. Here, under such condi- 
tions, we succeeded at phase connection with 8.4 GHz, but 
not with 2.3 GHz data. 

We removed most of the structure phase from the as- 
trometric observables based on the hybrid maps obtained 
from our observations and the selection of reference points 
in the maps, as in the past. Nevertheless, errors at the mil- 
liarcsecond level in the comparison of astrometric positions 
determined from data obtained at different epochs can re- 



sult from changes in the structures of the radio sources, if 
such changes are not taken into account properly. We try 
to diminish such error or bias, by attempting to remove the 
effects of the source structure from our data through refer- 
ring all astrometric measurements to the core (i.e., the same 
"physical" point) in the source. 

As a possible improvement to the general astrometric 
technique, we removed most of the effect of the ionosphere 
by means of TEC estimates deduced from GPS measure- 
ments at stations in the vicinity of the radio telescopes that 
we used. This important step forward is vital when ade- 
quate dual-band VLBI observations are not available. More- 
over, given sufficiently accurate corrections for ionospheric 
effects from GPS measurements, only single-band VLBI ob- 
servations need to be performed so that higher sensitivity 
(e.g., for fainter sources) can be obtained by using the whole 
recording bandwidth for this single band. The technical as- 
pects of the procedure will be presented in a complementary 
paper (Ros et al. 1999). 

After obtaining the pairwise relative positions of the 
three radio sources, we compared the separation of the 
pair QSO 1928-t-738/BL 2007+777 with two previous results 
from earlier epochs. Our relative position fell between the 
other two: we interpreted the differences as being due to dif- 
ferent "physical" reference points in the sources being used 
(unintentionally) at the different epochs of observation. 

According to H92, QSO 1928+738 is a superluminal ra- 
dio source that ejects components which seem to move along 
a sinusoidal trajectory starting at a flat-spectrum compo- 
nent detected at 22 GHz that H92 identified with the core. 
However, the inferences drawn from what seemed to be the 
correct alignment of our maps, led us to question this iden- 
tification, suggesting instead that the core is further north. 
We then identified four components that emerged during 
the six years spanned by our observations (late 1985 to late 
1991), which we labeled AO, Al, A2, and A3; their proper 
motions all appear to be in the range of 0.3-0.4 mas/yr. 

In the future, the improvements incorporated into our 
present data reduction can perhaps be successfully applied 
to radio sources much further apart on the sky and at any 
part of the solar activity cycle with the use of only single- 
band observations. In addition, the better performance in 
recent years of the global VLBI network (especially the 
VLBA) provides sizable advantages for phase-delay astrom- 
etry: higher slew speeds of radio telescopes with correspond- 
ing reduction of slew times, better sensitivity of receivers, 
and automated data correlation with the corresponding re- 
duction of "waiting" times. We believe that far more ambi- 
tious projects are now possible. 
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